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ABSTRACT 



As an extension of our previous work, we investigate the dynamical instability 
against nonaxisymmetric bar-mode deformations of differentially rotating stars in 
\^ ' Newtonian gravity varying the equations of state and velocity profiles. We performed 

the numerical simulation and the foUowup linear stability analysis adopting polytropic 
equations of state with the polytropic indices n = 1, 3/2, and 5/2 and with two types 
of angular velocity profiles (the so-called j-constant-like and Kepler-like laws). It is 
confirmed that rotating stars of a high degree of differential rotation are dynamically 
unstable against the bar-mode deformation, even for the ratio of the kinetic energy 
to the gravitational potential energy (3 of order 0.01. The criterion for onset of the 
. bar-mode dynamical instability depends weakly on the polytropic index n and the 

' angular velocity profile as long as the degree of differential rotation is high. Gravi- 

\ tational waves from the final nonaxisymmetric quasi-stationary states are calculated 

. in the quadrupole formula. For proto-neutron stars of mass 1.4M0, radius 30 km 

and /3 <^ 0.1, such gravitational waves have the frequency of ^ 600-1,400 Hz, and the 
effective amplitude is larger than 10~^^ at a distance of about 100 Mpc irrespective 
of n and the angular velocity profile. 

o 
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1 INTRODUCTION a central body with constant specific angular momentum 

are unstable against the Papaloizou-Pringle instability (Pa- 

: In our previous paper, we studied dynamical bar-mode sta- paloizou & Pringle 1984). On the other hand, the accretion 

- - ' bilities of differentially rotating stars in Newtonian gravity disks are stable if the velocity profile is Kepler-like. One 

(Shibata et al. 2002). In that study, we adopted a polytropic could claim that the bar- mode instability which we found 

equation of state with the polytropic index n = 1 and the would not be universal and might set in only for the very 

so-called "j-constant-like" angular velocity profile in which special rotational profile such as the j-constant-like law as 

the magnitude of the angular velocity decreases as at in the Papaloizou-Pringle instability. In addition, we focused 

large values of m, where ro denotes the cylindrical radius. We only on a stiff equation of state with n = 1 in the previous 

found that rotating stars of a high degree of differential ro- paper, so that one could also ask if the instability sets in for 

tation are dynamically unstable even for (3 = \T/W\ ~ 0.03, softer equations of state. 

where T and W are rotational and gravitational potential answer these questions, we have performed numer- 

energies. This value is much smaller than the long-believed j^al simulations of differentially rotating stars varying the 

criterion of ^ 0.27 for onset of the bar-mode dynami- polytropic index and the angular velocity profile. In this pa- 

cal instability of rotating stars (see Shibata et al. 2002 for pgr^ we report the numerical results. We will show that irre- 

a review). We also found that after the instability sets in, spective of the polytropic index n and the angular velocity 

such unstable rotating stars with 0.03 ^ /3 ^ 0. 15 eventually profile, a rotating star of a high degree of differential rotation 

settle down to nonaxisymmetric ellipsoidal quasi-stationary jg dynamically unstable against the bar-mode deformation 

states. even for the case that P is of order 0.01. We will also show 

However, there are two questions which have not been that an unstable star of a small value of /3 eventually settles 

answered in the previous work (Shibata et al. 2002). One down to a nonaxisymmetric quasi-stationary state, which is 

is associated with our choice of the j-constant-like angular a strong emitter of quasi-periodic gravitational waves, 

velocity profile. It is well known that accretion disks around The paper is organized as follows. In Section 2, we de- 
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scribe our methods in numerical analysis. In Section 3, the 
numerical results are presented. Section 4 is devoted to a 
summary and discussion. Throughout this paper, we use the 
geometrical units of G = c = 1 where G and c denote the 
gravitational constant and the light velocity. 



2 METHOD 

2.1 DifFerentially rotating axisymmetric stars 

We set a diflFerentially rotating star in equilibrium, and in- 
vestigate the dynamical stability. Rotating stars in equilib- 
rium are modeled using the polytropic equations of state 
as P = Kp^ where P, p, K and F = 1 + 1/n denote the 
pressure, density, polytropic constant and adiabatic index. 
In this paper, we choose n = 1, 3/2, and 5/2 (F = 2, 5/3, 
and 7/5). 

As the angular velocity profile 0,{w), we choose the so- 
called j-constant-like law as 



^2 + yl2 ' 



and the so-called Kepler-like law as 



O = no 



3/4 



(1) 



(2) 



where ^ is a constant, and ilo the angular velocity at the 
symmetric axis. The parameter A controls the steepness of 
the angular velocity profile: For smaller values of A, the pro- 
file is steeper and for ^ — > oo, the rigid rotation is recovered. 
In the present work, the values of A are chosen among the 
range Q.l < A = A/Rcq < 1 whore _Roq is the equatorial 
radius of rotating stars. For the rotation laws (1) and (2), 
Q at large cylindrical radius asymptotically behaves as 
and zu^^^^. This is the reason why we refer to the profiles 
(1) and (2) as the ji-constant-like and Kepler-like laws. 

In the limit of ^ ^ with the profile (1), the specific 
angular momentum becomes constant everywhere and f2 di- 
verges at 1*7 = 0. We note that this profile has been often 
used in studies of nonaxisymmetric instabilities in tori and 
annuli (Papaloizou & Pringle 1984; Goodman & Narayan 
1988; TohUne & Hachisu 1990; AndaUb & TohUne 1997). In 
this paper, however, we do not consider tori and annuli and 
focus only on spheroidal stars for which the density is not 
zero at m = 0. Thus, wo cannot adopt this limiting profile. 

In terms of /3 = r/|VK| and A, one rotating star is de- 
termined for a given rotational profile and polytropic index. 
Thus, in the following, we often refer to these two parame- 
ters to specify a rotating star. Here, T and W are defined 
as 



W=^ f d'xp^, 



(3) 



(4) 



where p and are the mass density and the Newtonian grav- 
itational potential. To specify a particular model, we may 
choose the axis ratio of the rotating stars Ca instead of /3. 
Here Ca is defined as the ratio of the polar radius Rp to the 
equatorial radius JZeq, i.e., Ca = Rp/Req- For the equations 
of state and the angular velocity profiles that we study in 



this paper, the value of Ca monotonically decreases with in- 
crease of P for a given set of A and n. This is the reason 
that Ca can be a substitute for /3. 

2.2 Dynamical stability investigation 

To investigate the dynamical stability against nonaxisym- 
metric bar-mode deformations, we have performed the nu- 
merical simulation as well as the linear stability analysis. 
Below we explain the methods for our numerical computa- 
tion separately. 

2.2.1 Numerical simulation 

In the hydrodynamic simulation, we initially superimpose a 

nonaxisymmetric density perturbation to an axisymmetric 
equilibrium star. We focus mainly on a fundamental bar- 
mode, and simply add the node-less density perturbation of 
the form 



5p = 6 ■ po{zu,z)- 



y 



(5) 



where po{vj, z) denotes the density of the axisymmetric con- 
figuration and 5 constant. Throughout this paper, we choose 
5 = 0.1. For simplicity, the velocity is left to be unperturbed 
at t = 0. The growth of the bar- mode can be followed by 
monitoring the distortion parameter as 



where 



Ixx ~\~ ^yy 

"■*-xy 
I XX H~ lyy 



(6) 

(7) 

(8) 



and lij {i,j = x,y,z) denotes the quadrupole moment de- 
fined by 



In = I d'^xpx'x^ . 



(9) 



Here, = {x,y,z). Simulations are performed using a 
3D numerical hydrodynamic implementation in Newtonian 
gravity (Shibata et al. 1997) (see also Shibata 2000 for re- 
sults of various test simulations with the identical hydrody- 
namic numerical scheme but in general relativity) . We adopt 
a fixed uniform grid with size 141 x 141 x 141 in x — y — z, 
which covers an equatorial radius by 50 grid points initially. 
We also performed test simulations with size 71 x 71 x 71 
(i.e., the grid spacing becomes twice larger) for several se- 
lected cases and confirmed that the results depend weakly 
on the grid resolution. Wo assume a reflection symmetry 
with respect to the equatorial plane. Since several rotating 
stars that we picked up have a flattened conflguration, we 
set the grid spacing of z half of that of x and y. 

2.2.2 Linear stability analysis 

In the linear stability analysis, wo employ the scheme de- 
veloped by Karino ct al. (2000, 2001). The Euler pertur- 
bations of the physical quantities are replaced by functions 
of the form /(r, 6»)e*""^-*'^* in the the linearized hydrody- 
namic equations. Here m is an azimuthal mode number. As 
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Figure 1. The dynamical bar-mode stability is shown in the 
plane of /3 and A for the j-constant-like angular velocity profile. 
The circles and crosses denote that the rotating stars are unstable 
and stable, respectively. The dashed curves denote the boundary 
which distinguishes the spheroidal stars from the toroidal stars 
(spheroidal stars are located below the curves). 

a result, the problem reduces to the eigenvalue problem for 
an eigenvalue u) and the corresponding eigenfunctions of the 
perturbed quantities. We assume an adiabatic relation be- 
tween the Euler perturbation of the pressure and that of the 
density. Since the fundamental mode of the m = 2 or bar- 
type oscillations is node-less, we have checked whether the 
obtained eigenfunctions satisfy that condition. In this paper, 
we have analyzed the stability of the oquihbrium configura- 
tions of n — 1 polytropes for two rotation laws (1) and (2) 
with several values of A. 



3 NUMERICAL RESULTS 

3.1 Dynamical Bar-mode Stability 

3.1.1 Results by numerical simulation 

The dynamical stability is studied for various combination 
of r, A and (3, and for two angular velocity profiles in the 
numerical simulation. In Figures 1 and 2, we summarize the 
results with regard to the dynamical bar-mode stability for 
the j-constant-like and Kepler-like angular velocity profiles, 
respectively. Here, the circles (crosses) denote that the stars 
of a given set of A and /3 are dynamically unstable (stable). 
We focus only on the spheroidal stars which are located 
below the dashed curves plotted in Figures 1 and 2. (If the 
value of /? is larger than that on this curve for a given value 
of A, the star is toroidal.) 

It is found that many rotating stars of a high degree of 
differential rotation as A = 0.1 and 0.3 are dynamically 
unstable even for (3 of order 0.01. In the case of the j- 
constant-like angular velocity profile, most of rotating stars 
with P J; 0.01 and A = 0.3 arc dynamical unstable both for 
r = 2 and 5/3. In the case of the Kepler-like angular velocity 
profile, the threshold of A for onset of the instability for a 
given set of /3 and F appears to be slightly smaller than that 



Figure 2. The same as Figure 1 but for the Kepler-like angular 
velocity profile. 



in the j-constant-like case. This is because the Kepler-like 
angular velocity profile is not as steep as the j-constant-like 
one for the same value of A. Indeed, for F = 2 and 5/3 
and for A = 0.3, there is a wide parameter space around 
(5 0.15 in which the stars are stable against the bar-mode 
deformation in the Kepler-like case. Such a wide parameter 
space is absent in the j-constant-like angular velocity profile 
for A — 0.3. However, for A — 0.1, the stars in a wide range 
of P are unstable oven in the Kepler-like angular velocity 
profile. Thus, rotating stars of a high degree of differential 
rotation with /3 of order 0. 01 are dynamically unstable even 
in the Kepler-like angular velocity profile. 

The threshold of A for onset of the dynamical instability 
for a given angular velocity profile also appears to depend 
on F: For a smaller value of F, the threshold of A is smaller. 
For example, compare the results with F = 5/3 and 7/5 
in Fig. 2. For F = 5/3, the stars oi A = 0.3 and 0.03 
<J 0.1 are unstable, but for F — 7/5, all the stars with 
P < 0.2 that we studied are stable for A = 0.3. The reason 
is that the rotating stars of smaller values of F have a more 
centrally condensed structure and as a result the effective 
degree of differential rotation of the central core for the stars 
of soft equations of state can become very high only for a 
sufficiently small value of A. 

The rotating stars of a high degree of differential ro- 
tation arc also dynamically unstable for a high value of 
P (J; 0.2) (see the results of ^4 = 0.6 and 1 for the j-constant- 
like law and of A = 0.3 for the Kepler-like law). The axial 
ratio Ca of these stars is very small (Ca ^ 0.3). In particu- 
lar, for the unstable stars of the Kepler-like angular velocity 
profile with ^ 0.2 and A = 0.3, Ca ~ 0.15. Thus, they 
have an almost toroidal shape. On the other hand, for the 
unstable stars with a small value of {f^ 0.10), Ca ^ 0.5 
and, hence, the shape is not significantly toroidal. 

An interesting feature is found for F = 2 and A = 0.6 
of the j-constant-likc angular velocity profile, and for F — 
5/3, 2 and A = 0.3 of the Kepler-like angular velocity profile. 
In these cases, the stability does not change monotonically 
with increase of 0: (1) stars of a sufficiently small value of 
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P 0.01) are stable (we have not carried out simulations 
for P < 0.01, but since the value of Ca with /3 ^ 0.01 is 
larger than 0.9 (i.e., the star is almost spherical), we assume 
that the stars with (3 ^ 0.01 are stable), (2) stars of 0.01 ^ 
/3 <; 0.1 arc unstable, (3) stars of 0.1 ^ /3 ^ 0.2 arc stable 
again, and (4) stars of a sufficiently large value of (3 (^ 0.2) 
are unstable again. 

It has been widely believed that the value of /3 is a 
good indicator to distinguish unstable stars from stable ones. 
However, the examples shown here illustrate that (3 is not a 
good indicator for determination of the dynamical stability 
of rotating stars of a high degree of differential rotation. 

The reason for the fact that the stability docs not 
change monotonically with /3 is not clear. The likely reason 
is that the mode associated with the stability for the high 
value of (3 is different from that for the small value. However, 
this interpretation cannot be proved at present because of 
the following reason. The property of the unstable modes 
for the high value of such as the perturbed density profile 
and the magnitude itself (j3 ^ 0.2) is essentially the same 
as that of the well-known bar-mode (i.e, the m = 2 toroidal 
mode, see Chandrasekhar 1969). Hence, we identify the un- 
stable modes of high value of /3 as the / mode. On the other 
hand, for the small value of f3, the perturbed density profile 
is also node-less, and furthermore, the real part of the eigen 
frequency is approximately proportional to {M/Rlq)^^"^ (see 
Figures 10 and 11). These are the properties that the / mode 
should have. This fact suggests that the unstable modes for 
the small value of (3 might be also the / mode. 

We suspect that there is something different between 
two unstable modes of high and low values of j3. However, 
to find the difference, it is necessary to precisely understand 
the definition of the / mode against the bar-deformation for 
differentially rotating stars. At present, it is not clear for 
us what properties, besides the node-less density profile and 
the cigen frequency of order of magnitude ~ {M / R^^y^'^ , 
characterize the / mode. Namely, we do not know what to 
prove. Precisely defining the / mode for differentially ro- 
tating stars is beyond scope of this paper, and is left as a 
problem in future. 

Before closing this section, we address the following 
point. As pointed out by CentrcUa ct al. (2001), the dy- 
namical instability for m = 1 mode may set in for rotating 
stars of a high degree of differential rotation. In this paper, 
we have not found the evidence that such instability sets 
in. This seems to be due to the fact that we adopted only 
the stiff equations of state. Centrella et al. (2001) studied 
the instability for differentially rotating stars of a very soft 
equation of state with n = 10/3. In such a soft equation of 
state, the instability associated with the m = 1 mode plays 
an important role. However, this is not likely the case in the 
stiff equations of state. Indeed, recently, Saijo et al. (2003) 
have pointed out that the instability of the m = 1 mode sets 
in only for differentially rotating stars of very soft equations 
of state (n ^ 2.5) and of a large value of [3 (^ 0.15). 

3.1.2 Results by linear analysis 

In Figure 3, wc display the numerical results of the linear 
stability analysis for the same equilibrium models as shown 
in Figures 1 and 2 for the n = 1 polytrope. These figures 
show that the results by the numerical simulation agree well 




I ■ ■ • • • • • 1 

0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



r=2 

0.15 - stable 
(b) Kepler 

_ 0.1 - » : ' "" i 

5 e o unstablfe 

0.05 - ' I I ..-'-"'^ 



I ' • • • 1 

0.1 0.2 0.3 0.4 0.5 0.6 



Figure 3. The results of the linear stability analysis are shown 
in the plane of (3 and A (a) for the ;?-constant-like and (b) for 
the Kepler-like angular velocity profiles. The circles and crosses 
denote that the rotating stars are unstable and stable, respec- 
tively. The solid curve in the upper panel denotes the threshold of 
the dynamical instability against the well-known bar-mode which 
was calculated in Karino and Eriguchi (submitted). The dotted 
curves denote the approximate threshold of the stability for the 
small value of (3. 

with those by the linear stability analysis. This fact con- 
firms the conclusion that rotating stars of a high degree of 
differential rotation are dynamically unstable even for the 
small value of 3 irrespective of the angular velocity profile. 
We note that the solid curve in the upper panel of Fig- 
ure 3 denotes the threshold of the dynamical stability (i.e., 
above this curve the star is dynamically unstable) against 
the well-known bar-mode (Chandrasekhar 1969) which was 
calculated and reported by Karino and Eriguchi (2003). 

3.2 Fate of unstable stars 

In this section, we focus on rotating stars of the Kepler-like 
angular velocity profile, since the results for the j-constant- 
like angular velocity profile with F = 2 arc shown in the 
previous paper (Shibata et al. 2002) and, moreover, we have 
found that the results for F = 5/3 show qualitatively iden- 
tical features. 

In Figures 4-7, we display time evolution of as a func- 
tion of Q.Qt for (r,i) = (2,0.3), (5/3, 0.3), (5/3, 0.1), and 
(7/5, 0.1). As shown here, the value of rj does not reach the 
value of order 1 but saturates at order 0.1. This implies that 
the growth of the instability saturates at a weakly nonlin- 
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Figure 4. Time evolution of r; as a function of Qgt for F = 2 and 
A = 0.3, and for the Kepler-like angular velocity profile. 
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Figure 6. The same as Figure 4 but for F = 5/3 and A = 0.1. 
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Figure 5. The same as Figure 4 but for F = 5/3 and A = 0.3. 
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Figure 7. The same as Figure 4 but for F = 7/5 and A = 0.1. 



ear stage. We note that for unstable stars with f3 ^ 0.2 and 
A = 0.3, the value of r] increases to ~ 1 irrespective of F, 
implying that a highly deformed star is produced. There are 
a number of numerical works in which such a highly nonax- 
isyrnrnctric structure is the outcome after onset of dynamical 
instability of differentially rotating stars with a high value of 
/? (e.g., Williams & Tohline 1987, 1988; Houser & Centrella 
1996). The outcome in the simulations for unstable stars of a 
high value of /3 found in the present numerical computation 
is qualitatively the same as that in previous papers. Thus, 
we do not discuss the results for such cases in the following. 

After the perturbation saturates, the amplitude of r) 
settles down toward a value of order 0.1. The approximate 
final values of rj are 



Vf < 



~ 0.1 for {r,A) = (2,0.3), 

^0.1 for (r,i) = (5/3,0.3), 

~ 0.2 for (r,i) = (5/3,0.1), 

~ 0.03 for (r,i) = (7/5,0.1). 



(10) 



Comparing the results of {r,A) = (5/3,0.3) and (5/3, 0.1) 
for the identical value of Ca, it is found that for the smaller 
value of ^, r]f is larger. This implies that a stronger degree 



of differential rotation makes the magnitude of the nonax- 
isymmetric deformation larger. 

It is also found that for the smaller value of F, the final 
value of r? is smaller. This is because the stars of smaller 
values of F have a more centrally-condensed structure and, 
hence, their efTective steepness of the differential rotation is 
smaller for smaller value of F for a given set of Ca and A. 

In Figures 8-9, wc display the snapshots of the density 
contour curves in the equatorial plane at selected time steps 
for {T,A,Ca) = (2,0.3,0.7) and (5/3, 0.3, 0.8). The value 
of 3 is about 0.071 and 0.045, respectively. In both cases, 
the nonaxisymmetric perturbation initially provided grows, 
changing the shape of the rotating stars to be ellipsoidal. 
However, the perturbation does not grow to the highly non- 
linear stage and, hence, neither a spiral arm nor a large bar 
is formed in contrast to the outcome in the simulations with 
a high value of /3 ^ 0.2. Instead, the slightly deformed el- 
lipsoid is the final outcome. Since the deformed ellipsoid is 
almost stationary, 77+ and r)x oscillate quasi-periodically in 
the late phase of the simulations. This result is qualitatively 
the same as for the stars of the j-constant-like angular ve- 
locity profile (Shibata et al. 2002). 

In Figures 10 and 11, we show the frequency fr of the 
oscillation of the ellipsoidal star in units of {M/Rg^)^^^ as 
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Figure 8. The density contour curves at selected time steps for 
r = 2, A = 0.3, and Ca = 0.7 and for the Kepler-like angular 
velocity profile. Here, Pq is 2n/Qo. The contour curves are drawn 
for p/pmax = 0.95, 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2, 0.1 0.01 
& 0.001, where pmax denotes the maximum density at each time 
slice. The dashed curves are plotted for p/pmax = 0.01 and 0.001. 
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Figure 9. The same as Figure 8, but for F = 5/3, A = 0.3 and 
Ca = 0.8. 
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Figure 10. fr of the ellipsoids formed after onset of dynamical 
instability for the j-constant-like angular velocity profile. Filled 
circles and open squares denote the results for (F, A) =(2, 0.3) 
and (5/3, 0.3). 
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Figure 11. The same as Figure 10, but for the Kepler-like angular 
velocity profile. Filled circles, open squares, filled triangles and 
crosses denote the results for (F,A) =(2, 0.3), (5/3, 0.3), (5/3, 
0.1) and (7/5, 0.1). 



a function of initial values of P for various sets of F and A. 
The value of fr is determined by the Fourier transform of ri+ 
in the time domain. We note that the rotational period of 
the ellipsoid is 2//,., and that the frequency of gravitational 
waves is fr. It is interesting to note that for F = 2 and 5/3, a 
nondimensional quantity fr = fr{R'cq/M)''^^ is in a narrow 
range between 0.2 and 0.35 irrespective of A, (3 and angular 
velocity profile. For F = 7/5, fr is between 0.4 and 0.55, 
which is larger than the values for F = 2 and 5/3. However, 
it is still in a narrow range. The fact that the value of fr 
depends weakly on (3 suggests that the excited mode may 
be the / mode. 



3.3 Gravitational waves 

As discussed in our previous paper (Shibata et al. 2002), the 
dynamically unstable rotating stars, which are deformed to 
nonaxisymmetric ellipsoidal objects, are likely to be sources 
of quasi-periodic gravitational waves. In Figures 12-14, we 
show the gravitational waveforms of the -I- mode along the 
z-axis (/i+) s^nd the luminosity (E) as a function of time for 
(F, A) = (2, 0.3), (5/3, 0.1) and (7/5, 0.1). For aU the models 
picked up here, Ca = 0.705. The value of /? is fs 0.071, 0.046, 
and 0.047, respectively. We note that the waveforms for the 
X mode are essentially the same as those for the -|- mode 
except for the phase difference by 7r/4. Here, we calculate 
gravitational waves in the quadrupole formula (Misner et al. 
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Figure 12. Gravitational waves for /i+ in units of Af^/Rcq and 
the luminosity of gravitational waves E in units of (M/RcqY' as a 
function of Qot for (T, A) = (2, 0.3) and Ca = 0.705 (f3 Ri 0.071). 



Figure 13. The same as Figure 12, but for {r,A) = (5/3,0.1) 
and Ca = 0.705 (/3 « 0.046). 



1973) and define the waveforms by 



2L 



r 



and the luminosity by 



(11) 



(12) 



where 



JL - T T 



d li. 



r(3) 



and r is the distance from a source to a detector. /i+ and 
/ix are the waveforms observed along the z axis. 

In the early phase before the growth of the nonax- 
isymmetric perturbation saturates, the amplitude of grav- 
itational waves increases gradually together with the mag- 
nitude of r\. Then, the growth of the amplitude saturates and 
subsequently, quasi-periodic gravitational waves are emitted 
for a much longer time than the rotational period. The am- 
plitude of quasi-periodic gravitational waves depends on F 
and 4, as in the case of the saturated value of r). For a 
given set of A and Ca (or /3), it is smaller for smaller values 
of F. This is because the star with smaller values of F is 
more centrally condensed and, hence, the magnitudes of the 
quadrupole moments are smaller for the identical mass and 
radius. 

Wc found that for the parameters we studied (/3 0.1, 
A < 0.6 and angular velocity profiles (1) and (2)), the 
typical values of the amplitude and luminosity of gravita- 
tional waves for F = 2 are rh+.x ^ 0.1 0.2(A/'^/_Roq) and 
E ~ 0.005-0.01(Af/i?eq)*. For F = 5/3, the magnitudes of 
them are only slightly smaller than those for F = 2. How- 
ever, for F = 7/5, they are smaller by a factor of several as 
r/i+,x ~ 0.05(A/V-Rcq) and E ~ 0.001(M/iicq)^ 

Using the results shown in Figures 12-14, wc can esti- 
mate the expected effective amplitude of gravitational waves 
from the nonaxisymmetric outcomes formed after the dy- 
namical instability saturates. Here, we pay particular at- 




400 600 

n„ t 



800 1000 



Figure 14. The same as Figure 12, but for (T,A) = (7/5,0.1) 
and Co = 0.705 (/3 » 0.047). 



tention to proto-neutron stars which are likely formed soon 
after the supernovae of mass ~ 1.4M0 and radius several 

10 km. For the luminosity E = 0.001e(Af/_Rcq)'^ where e is 
a parameter of magnitude 1-10, the emission timescale of 
gravitational waves can be estimated as 



r^l = 100aoe-(A 



= 6.1 sec CKoe 



0.1 



M 



30 km 



M 



I AM, 



-3 

14) 



where we set T — aof3M^ /Req, and ao is a constant which 
depends on the value of F but very weakly on A; for /3 ;^ 0.1, 
ao ~ 0.8, 0.9, and 1.2 for F = 2, 5/3, and 7/5, respectively, 
within ~ 10% error. 

The characteristic frequency of gravitational waves is 
denoted as 



/ = /. 



: 790 Hz (A 



30 km 



-3/2 



M 



I AM, 



1/2 



(15) 



Assuming that the nonaxisymmetric perturbation would not 
be dissipated by viscosities or magnetic fields on the emis- 
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sion timescale of gravitational waves (e.g., Baumgarte et al. 
2000), the accumulated cycles of gravitational wave-train N 
are estimated as 



N = fT 



Re, 



0.1 / V 30 km 



5/2 



M 



lAMr: 



-5/2 



(16) 



The effective amplitude of gravitational waves is defined by 
hes = N^^^h where h denotes the characteristic amplitude 
of periodic gravitational waves. Using this relation, we find 



30 km 



1/4 



M 



I AM, 



3/4 



100 Mpc 



(17) 



(Thorne 1987; Lai & Shapiro 1995; Liu & L_indblom 2001; 
Liu 2002) where h = hrReq/M^. Since fr, h, e and /3 de- 
pend on the values of F, A, and Ca, h^s can vary by a factor 
of --^ 3. However, for all the rotating stars that we studied 
in this paper, hes is always larger than 10^^^ at a distance 
r ~ 100 Mpc with i?eq ~ 30 km and M ~ IAMq. Fur- 
thermore, the frequency of gravitational waves is about 1 
kHz for _Roq ~ 30 km and M ~ IAMq. Thus, gravitational 
waves from proto-neutron stars of a high degree of differen- 
tial rotation, of mass ^ IAMq, and of radius ^ 30 km at a 
distance of ~ 100 Mpc are likely to be sources for laser in- 
terferometric detectors such as LIGO (Thorne 1995), if the 
other dissipation processes are negligible. 



4 SUMMARY AND DISCUSSION 

Wo have studied the dynamical bar-mode instability of dif- 
ferentially rotating stars of polytropic equations of state. 
We chose three polytropic indices and two angular veloc- 
ity profiles in this study. We found that rotating stars of 
a high degree of differential rotation are dynamically un- 
stable against the nonaxisymmotric bar-mode deformation 
even with (3 <C 0.27, irrespective of the polytropic indices 
and angular velocity profile. The criterion of the value of (3 
for onset of the instability depends on the rotational pro- 
file and the equations of state, but the dependence is very 
weak if the degree of differential rotation is high enough as 
A - 0.1. 

We estimated the effective amplitude of gravitational 
waves from nonaxisymmotric objects formed after onset of 
the dynamical instability. For typical proto-ncutrori stars of 
mass ~ IAMq and radius several 10 km, the cfi^cctive am- 
plitude of gravitational waves at a distance of ~ 100 Mpc is 
larger than 10~^^ , and the frequency ~ 1 kHz. Therefore, the 
gravitational waves can be sources for laser interferometric 
detectors such as advanced LIGO (e.g., Thorne 1995). 

As we mentioned above, this conclusion is drawn under 
the assumption that dissipation of nonaxisymmetric pertur- 
bations by viscosity and magnetic fields is negligible. The 
dissipation timescale due to molecular viscosities and mag- 
netic braking is likely to be longer than 10 sec (e.g., Baum- 
garte et al. 2000). Thus, these effects can be safely neglected. 
However, turbulent magnetic viscosity (Balbus & Hawley 



1998) may be relevant for redistribution of the angular mo- 
mentum profile. This implies that a differentially rotating 
star might be enforced to a rigidly rotating state on a dy- 
namical timescale and, hence, the nonaxisymmetric struc- 
ture might disappear. Many theoretical works have clarified 
that the magnetic viscous effect can redistribute the angu- 
lar momentum distribution of accretion disks around central 
objects on a dynamical timescale (Balbus & Hawley 1998). 
However, to our knowledge, there is no work on the magnetic 
field effect to sclf-gravitating rotating stars. The magnetic 
field would redistribute the angular momentum of the differ- 
entially rotating stars on a dynamical timescale, but it is not 
clear if it is strong enough to enforce the rotating stars to 
an axisymmetric state in a dynamical timescale. To clarify 
this problem, it is necessary to carry out MHD simulations 
for rotating stars. Such work should be done in the future. 
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